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Abstract 



We compute, in the MSSM framework, the total electroweak contributions at one loop for the 
process pp — > tW + X, initiated by the parton process bg — * tW. The supersymmetric effect 
is analyzed for various choices of the SUSY benchmark points. Choosing realistic unpolarized 
and polarized experimental quantities, we show the size of the various effects and discuss their 
dependence on the MSSM parameters. 
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I. INTRODUCTION 



The relevance of the process of associated tW production from proton - proton collisions 
has been exhaustively stressed in recent dedicated studies Q|. In the Standard Model frame- 
work, it is well known that accurate measurements of the production rate would provide an 
excellent determination of the tbW coupling. For physics beyond the SM, one expects that 
precision tests of virtual effects might be performed, provided that the effects were suf- 
ficiently large, i.e. at least of the same size as the overall theoretical and experimental 
uncertainty. On the theoretical side, an estimate given in Q| predicts for the total produc- 
tion cross section an uncertainty of about 15%. On the experimental side one must recall 
the fact that the considered process will be seen, for the first time, at LHC, simply because 
of the required pp energy. Therefore an estimate of the expected experimental uncertainty 
is in fact still missing. This might be particularly relevant, if the estimated effect turned 
out to be reasonable (e.g. of the same size as the theoretical one), for the special purpose of 
performing a precision test of supersymmetric models, in particular of the simplest available 
one, the MSSM. In fact, in a previous paper jsj, the genuinely weak effects of the model 
were considered at one loop for all the processes of single top production (td, tb, tW and also 
tH~) at LHC. The treatment was rather preliminary and qualitative, and only considered 
the very special case of a light SUSY scenario and of a production in a large (~ 1 TeV) 
final invariant mass range, where a simple logarithmic expansion of so called Sudakov kind 
could be used. The feature that emerged was that, for the three processes that will certainly 
be seen at LHC (i.e. td, tb, tW), the electroweak relative effect in the MSSM is sizable, 
particularly for large tan/3 values where it could reach the 30 % size. This appeared to us a 
good motivation for performing a complete accurate estimate, valid for all realistic invariant 
masses and containing all the parameters of the model. In this paper, we present the results 
of our effort for the tW production process. It is, to our knowledge, the first complete cal- 
culation of the electroweak MSSM effects at one loop for the process, that also includes the 
QED soft photon radiation. We have checked the validity of our results in three different 
ways, i.e. (a) we have verified the cancellation of all the (virtual) ultraviolet divergences, 
(b) we have verified the cancellation of all the (real and virtual) infrared divergences and, 
last but not least, (c) we have verified the exact reproduction of the asymptotic Sudakov 
expansions, given in |2|, from the computed Feynman diagrams. After these three checks 
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we hope that our results should be correct, and we shall show them in the paper with the 
following plan: in Section |H] a brief description of the relevant Feynman diagrams is given, 
and a discussion of the cancellation of the ultraviolet and infrared divergences is provided. 
Section ITTT1 shows the reproduction of the (essentially academic) asymptotic Sudakov expan- 
sion. In Section II VI the realistic observables are defined, and the related MSSM effect is 
shown for various choices of the SUSY benchmark points. A final discussion that included 
a review of future calculations is provided in Section V. 



II. MSSM bg -► tW PRODUCTION AT ONE LOOP 

The process that we have considered is the so called exclusive associated tW production, 
whose partonic description corresponds to the two body final state reaction 

bg - tW~ (1) 

that is represented, at Born level, in Fig. (1). In this paper we shall not consider the inclusive 
process 

gg^tbW-. (2) 

As known the collinear b component of this process is already enclosed as a QCD NLO 
correction to the bottom quark distribution function of the exclusive process Eq. (1) and 
our treatment will only consider the one-loop electroweak effects. 

At the one loop level, we have to consider different kinds of Feynman diagrams, several 
of which will exhibit an ultraviolet divergence. We shall choose the on-shell renormalization 
scheme, and in this framework we shall define the following classes of Feynman diagrams: 



A. Born, self-energies and counterterms 

The two Born diagrams represented in Fig. 1 (a, b) are an s-channel b quark exchange 
and a w-channel top quark exchange. With the definitions s = (p& + p g ) 2 = {pw + Pt) 2 , and 
u — (Pb — Pw) 2 = (Pg — Pt) 2 channel, we have 



.JP L { g + m b )t , M'+m t )^P L 

where 



A»™(gb - wt) = -^a(t) r yn 2 + — 7 Mb) (3) 
s w V2 s-m 2 b u-m 2 t J 



q = p 9 +Pb=Pw+Pt s = q q' = p t - p g = p b - p w u = q , (4) 



and e, e are the gluon and W polarization vectors, respectively. 

In the on-shell renormalization scheme, these Born terms have to be completed with 
counterterms associated to the b, t, and W lines. These counterterms are expressed in 
terms of quark and gauge bosons self energy functions Tj^ LRS (k 2 ), S yv "(/c 2 ). In these self- 
energies we take into account SM and SUSY contributions (sfermions, Higgs, neutralinos 
and charginos). The b and t quark propagators are also modified by self-energy functions 
of s and u. 

In the s-channel, we can use the invariant forms 



and write the amplitude as 



where 



Iil.r = WPl,r IL,r = &Pl,r (5) 



E W % + n° "/y (6) 

v 



n? L = — -{i + 5zY - sz? + hv w + h* t 

\/2sw{s - mi) 2 2 

+pzl + \sz' L - j^VM + szi) 



N{ R = iV 2 s L = (8) 



iv 2 s R = j= eg ; mb 2 M + 5zY - szY + \w w + h* t 

\J2s w {s - mi) 2 2 



+ ^-5Z b L + ^5Zl + 5Z b R sQ* a (s) + 5Z b R + E* (*) + SZ b L ) 

2 2 s — mf 



1 „ „>. 1 , „ „fc S 

In the w-channel, we define 

I\L,R = MtPL,R l2L,R = tiPL,R (10) 

and write 

A = E W "i?,, + N u 2 V" } (11) 
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with 



L = -=-^1 -{1 + 5ZY - 6Z? + \s* w + ±M t 

V2s w {u - mf) 2 2 

+hzl + \8Zl - -^— 2 {Y* L {u) + 5Z{) 

Zi Zi LL /I l-f^ 

m ^n(n) + SZfo - r^PsM " + Mr) ~ ^1) (12) 



u-m 2 t y H/ u-m\ v * K ' 2 y ^ u/ m t 



N? R = N'2 R = (13) 

NX L = -= e9 ; m - {1 + 8zy - + \s* w + 1-5% 

V2sw{u — mf) 2 I 

1 1 11 

+-SZ b L + -5Z{ + 6Z* R - ——-^(u) + 5Z R + E* L (u) + 



i2 



U .- m? '^("»-2(^ + ^)-^fl> < M » 

The various counterterms have the following explicit expressions in terms of self-energies. 
First, we have the divergent quark wave function renormalizations 

5Z b L = SZl = 5Z L = -ElK 2 ) - ml^'Kml) + E^(m 2 ) + 2E' 5 V fe 2 )] (15) 



6Z b R = -E^(m 2 ) - m 2 [E2K) + E^(m 2 ) + 2E' 5 6 (m 2 )] (16) 

5Z R = 5Z L + T, t L {m])-T, R {ml) (17) 

Then, we have the finite wave-function renormalization required in the on-shell scheme and 
unavoidable since we have both up and down type quarks in the process 

6% = -{ElK) + 5Z L + m 2 [E2K) + E^rn 2 ) + 2E'j(m 2 )]} (18) 

The similar terms for the W gauge boson are 



sz^-sz™ = - v ;; 2 (19) 



E^ z (0) 



^ = _^ (0 ) + 2^'(0) + J[^-^l (20) 
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and 

6* w = -T! ww (M^) = -{£ W (M^) + 5Zf) (21) 
Finally, we list the mass counterterms 

6M^ = ReZ ww (M^) 5M 2 Z = ReZ zz (M 2 z ) (22) 

5m b = ^Re[E b L (ml) + E^(m 2 ) + 24(m 2 )] (23) 

5m t = ^ifepKm?) + S^K 2 ) + 2£< s (m 2 )] (24) 
B. Vertex corrections and Box diagrams 

The next two classes of diagrams are triangle-like vertices and box diagrams. A list of the 
generic diagrams (i.e. diagrams with virtual particles left unspecified apart from their spin) 
is shown in Fig. (J2J) as produced by FeynArts |4j. Schematically we can further subdivide 
them as follows (q stands for b or t quarks, V for 7, Z, W, H for neutral or charged Higgses 
or Goldstone particles, x f° r chargino or neutralino): 

1. Initial s-channel triangles connected to the intermediate b quark: (Vqq), (Hqq), (xQQ)] 

2. Final s-channel triangles connected to the intermediate b quark: (btV), (HHq), (btx) 
(VVq), (HVq), (VHq), (btH), 

3. Up it-channel triangles connected to the intermediate t quark: (qqV), (qqH), (qqx)i 

4. Down w-channel triangles connected to the intermediate t quark: (tbV), (tbH), (xXQ)> 
(VVq), (VHq), (HVq), (HHq), (tb X ); 

5. Direct boxes: (bbtx°), (bbtV), (bbtH); 

6. Crossed boxes: (qqVV), (qqVH), (qqHV), (qqHH), (qqxx)', 

7. Twisted boxes: (ttVb), (ttHb), (ttx°b). 
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The notation corresponds to the clockwise ordering of the internal particles inside the 
diagrams. 

An essential step consists then in checking the cancellation of UV divergences. They 
appear in the self energy functions S(fc 2 ), in the various counterterms and in the various tri- 
angles. Box contributions are convergent. We have checked the cancellation when summing 
all of these terms. This cancellation occurs in several independent sectors (gauge, Higgs, 
SM, SUSY). 

Having completed the first important check (cancellation of UV divergences) we now 
move to the forthcoming issue of cancellation of IR divergences that will be treated in the 
forthcoming discussion. 

C. Cancellation of IR divergences 

QED radiation effects are usually split into a soft part containing the potential IR singular 
terms, and a hard part including the emission of photons with energy not small compared 
to the process energy scale. In this brief section, we discuss the soft emission and the 
detailed cancellation of IR divergences that occurs when it is combined with virtual photon 
exchanges. 

Let us denote by A Born and A lloop any invariant helicity scattering amplitude evaluated 
at Born or one loop level. Let us also denote by A the photon mass acting as an IR regulator. 
The IR cancellation between (soft) real radiation and virtual photon exchange holds in every 
helicity channel separately and we have checked it numerically. 

It reads 

(A Born ) 2 (l + ^-5 S ) + 2A Born A lloop = finite as A -> (25) 

where, in the above expressions, 5s is the correction factor taking into account the emission 
of soft real photons with energy from A up to E™ ax « y/s. The explicit expression for 5,s 
can be found, for instance, in [5J]. 

In practice the above relation follows from the eikonal factorization 

A lloo P = _ A BornU §s + regular terms as _\ _> Q (26) 

4:71 

It is possible to split further the above factorization property. Indeed, the singular part of 
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the radiation factor has the form 

Ss = log £^ + re g ular terms as A -> (27) 

where i and j runs over the initial/final charged particles, i.e. (bt), (bW), and (tW). There 
are two types of contributions: the diagonal ones with i = j and the off diagonal ones with 

Now, the matching between the singular log A in the l.h.s. and r.h.s. of Eq. (j2*Uj) can be 
checked in several independent steps as follows 

1. the diagonal radiation terms i — j match the IR divergence in the counterterms 
associated to the i-th external line p. 

2. the off-diagonal radiation terms % ^ j match the IR divergence in the diagrams which 
are obtained connecting in all ways the i-th and j-th external lines with a virtual 
photon. This operation produces both triangle and box diagrams. 

As a final comment, we remark that gauge invariance is crucial to cancel all non factoring 
contributions associated to the final W line as discussed in |zj]. 

The next step in the treatment of QED effect is the calculation of hard photon emission. 
We have left this subject to a dedicated study which shall be discussed separately 



III. SUDAKOV EXPANSION OF THE SCATTERING AMPLITUDES 

Let us now consider the high energy behavior of the bg — > tW helicity amplitudes Fx^y^, 
where A, /i, A', // refer to the helicities A&, X g , A t , \w respectively. Several simplifications ap- 
pear in the Born and in the one-loop contributions. When s 3> m\ (m, being the internal or 
external involved masses), ignoring mf/s contributions, the non-suppressed Born amplitudes 
reduce to F , for transverse W and -F-++o for longitudinal W. 

The leading high energy Born helicity amplitudes are 



^Born 



^74>A (28) 
s w v 2 z cos | 



S 



pBom eg, \ l m t 9 
-+ +0 s w l 2 j M w 2 l l + cos0 j (oU) 



Note that -F-++o is controlled by the top Yukawa coupling factor ~ m t /Mw- In fact the 

amplitude o also occurs but at a much weaker level as it is controlled by the bottom 

Yukawa coupling factor ~ rrib/Mw- 

At one loop, these amplitudes receive logarithmic enhancements as discussed in several 
papers, called Sudakov terms. These terms are separated into universal and into angular de- 



pendent components. From the ru 
(there are misprints in the paper 
For transverse W amplitudes: 



es established in |9j, one expects the following expressions 
The correct equations are the following Eqs. (31-46)). 



1US 

2 ]: 



FY™,» = F*°™J - ( c™(bb) L + c ew (tt) L ) + c™(W T )] (31) 



c ew (qq)L = c ew {qq) L = c(qq, gauge) L + c(qq, yuk) L (32) 



c(dd, gauge) L = c(uu, gauge) L = a j]. + <2 f°^ > (n log— ^ log 2 — ^-) (33) 

— Ot S n S 

c(dd, gauge) R = ——^ (n log-^- - log —^) (34) 

c(uu, gauge) R = - — — (n log — r - log — r ) (35) 
where n = 3, 2 in SM and MSSM, respectively. 

2 2 

— — (X S Tfl Tflt 

c(bb, yuk) L = c(ti, yuk) L = - — — - [log — r ] [-fy t + -^y h \ (36) 

2 

c{bb, yuk) R = - [log -^-] [^-y b ] (37) 

2 

— OL S Tfl 

c(ti, yuk) R = - [log -J-] [-f y t ] (38) 



9 



where y t = 1, 2(1 + cot 2 /?) and y b = 1, 2(1 + tan 2 (3) in SM and MSSM, respectively. 



Oi r , 2 S 



and for the longitudinal W Q amplitude: 



<? + ,o = ^-,?, n + ,o[ \ ( ^(&% + c™(tt) R ) + d-(Wb)] (40) 



with, in SM: 



such that 



+ (41) 



Ct , r r , 9 S lr 13 + 14C? 



l + 2cgy m|_ 

[ 2s 2 r 2 2s 2 r 2 JL 8 M 2 J 1 1 ' 



whereas in MSSM : 



4tt u Ss^c 2 , b M^ J 

+P°g h + 7^(1 + tan 2 + -j^U + cot 2 /?)]} 



M 2 / JL 36*^ 4^71^ ^ 4^71^ 



(43) 



such that 



^-^I^KI-log'^ll^]} (44) 



(in which all single logs cancel !). 

For the electroweak angular terms we find: 



FT^ = F_T-,J- ^Hlog^K [log^liy^] + ^log^ } (45) 

z,n 1V1 W b 0{ -> b w c W b W b 
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F- + , = F-, [- ^lHog^K [flog^] - ^log^ } (46) 

Note that the longitudinal amplitudes satisfy the equivalence theorem which states 
that, neglecting mf/ s contributions, they should coincide with the amplitudes for the process 
bg — > tG~ , G~ being the charged Goldstone boson. 

We have checked, by using the asymptotic expansions of the B, C, D functions appearing 
in the self-energies, triangle and box amplitudes that our full one-loop result produces the 
logarithmic contributions expected by the rules given above. 

These resulting asymptotic expressions deserve several comments. In the case of trans- 
verse W production, one checks that at Born and one-loop level and at next-to-leading 
logarithmic accuracy in addition to trivial fermion chirality constraint \ t = A& = —1/2 
gauge boson helicity conservation [10] is preserved, both in SM and MSSM cases, i.e. only 
H = fi' amplitudes survive. One then sees that the MSSM differs from the SM in the sin- 
gle logarithm contributions, n = 2 instead of n = 3 for gauge terms and 2(1 + cot 2 /?) or 
2(1 + tan 2 /?) Yukawa enhancements, especially large for large tan/?. 

In the case of longitudinal W production, the Born amplitude is controlled by the Yukawa 
m t /My/ factor associated to fermion chirality violation A t = — A& = 1/2 and satisfies also 
the rule X g + A& = Xt which is an extension of the GBHC rule [10]. An additional remarkable 
feature appears for the single log contribution, namely it totally cancels in the MSSM case. 

Having successfully performed the ultraviolet, infrared and Sudakov tests, we hope that 
our complete expressions will be correct. In this respect, we should add the following 
comment: We do not expect that, at lower energies and for higher SUSY masses, the simple 
features that we met in the light SUSY Sudakov description given in |2| retain their full 
validity. Still, we would expect that, at least, some of the main features could survive. For 
instance, for what concerns the slope of the invariant mass distribution, we could hope that 
a simple modification at lower energies, or at lower energy/SUSY masses ratios, might be 
the addition of a (possibly large) constant term at least in a moderate energy region not 
too far from the asymptotic one, so that a smooth connection between the two regions is 
achieved. In the following section we shall return on this point, but first we shall define and 
examine those quantities that will be the realistic experimental observables. 
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IV. PHYSICAL PREDICTIONS 



We are now able to provide numerical predictions for the complete electroweak effect of 
the MSSM at one loop on the realistic observables of the considered tW production process. 
With this aim, we shall divide our presentation in two parts, that correspond respectively to 
the consideration of unpolarized and of polarized quantities. Following a pragmatic attitude 
i.e. assuming that only unpolarized observables will be measured in a first stage of the 
experiments, we shall start our analysis with the former ones. 

The first quantity that we shall consider is the invariant mass distribution, conventionally 
defined as 

da(PP ^ tW~ + X) 1 f c ° s9 ™* da 



/•cos V max J 

/ d cos 9 [L bg (r, cos 9) f^- {s)) (47) 

•/cos0 mi „ acostf 



ds S 
where r = ■% , and L bg is the parton process luminosity. 

L bg (r, cos 9)= dy [ b(x)g(-) + g(x)b(-) } (48) 

" ymin 

where S is the total pp cm. energy, and i(x) the distributions of the parton i inside the 
proton with a momentum fraction, x = related to the rapidity y of the tY system jll|. 

The parton distribution functions are the latest NNLO MRST (Martin, Roberts, Stirling, 
Thorne) set available on 3]. The limits of integrations for y depends on the cuts. We have 
chosen a maximal rapidity Y = 2 and a minimum pt which we shall specify later. 

Note that we are at this stage considering as kinematical observable the initial partons 
cm. energy i/i, and not the realistic final state invariant mass M t w- The transition from 
the first quantity to the second one can be performed using the available suitable event 
generators, 
production 

transition from yfs to M t w- This correction can be considered as a QCD effect, and as such 
it will be consistently treated in a forthcoming paper [jjj where this type of non electroweak 
effects will be included. For what concerns the complete one-loop electroweak amplitude, we 
can compute it for any choice of the MSSM parameters, but before doing this we want to show 
some features of the simple Born approximation of the partonic amplitude that we consider 
particularly relevant for an understanding of our following results. More precisely, the point 
that we want to stress is that the partonic invariant scattering amplitude for the process, that 
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ike for instance PYTHIA Q], as we did in a previous paper on top-antitop 
~\\ . We expect from that experience a small (few percent) modification in the 



represents the starting block of our calculations, turns out to be the sum of twelve different 
helicity amplitudes, that have been defined already in Section IIHI For large values of \/s, 
i.e. for ,/s sufficiently larger than the masses of all the particles and sparticles involved 
in the one-loop description of the process, we expect that only three helicity amplitudes 

remain dominant, more precisely those that have been defined in Section UTTl as F , 

-F_++o (the third and fourth index specifies the top and W helicity). The remaining 
helicity amplitudes vanish asymptotically i.e. for s — > oo like 1/s with possible logarithmic 
corrections at one loop, and in our preliminary paper |2] they were systematically neglected 
in the region of that was considered, corresponding to a) energies in the 1 TeV range and 
b) light SUSY masses scenario. For the realistic analysis that we can now carry on, both 
assumptions will be abandoned. In particular: 

a) The possibility of identifying the final (t, W) signal must face the serious competition 
of a background, mostly due to events coming from the copious top-antitop and WWj pro- 
duction. This problem has been already exhaustively discussed in a previous paper jjj, where 
it has been shown that the introduction of suitable b-tagging cuts will allow to extract the 
signal at reasonable (20 fb _1 ) luminosities. A priori, one would expect that the background 
contamination should be under control for cm. energies below a qualitative background 
threshold of approximately, say, 400-500 GeV, and increase in the higher energies region. 
Keeping this limitation in mind, we have nonetheless analyzed in this paper the full energy 
region from threshold to 1 TeV, although at this final energy value the identification of the 
signal might be difficult. The reason of this (optimistic) choice is that we do not have yet 
at disposal a rigorous experimental analysis of the realistically expected size of the signal 
at variable energies, as we had in the preliminary top-antitop paper 14j. This analysis is 
being already performed, and will be included in the already mentioned forthcoming work. 

b) The SUSY scenario that we shall investigate is the conventional mSUGRA one. In 
particular, we shall consider a number of benchmark points that are nowadays available, 
trying to choose those that show a definite difference in the values of the various SUSY 
masses, and of tan/3. We insist on the fact that we could perform its calculation for any 
choice of the parameters, but for obvious reasons we have limited the presentation of Figures 
in this paper. 

After these preliminary remarks, we now show in Fig. (J3J) the comparison (treated in Born 
approximation) of quantities that we consider particularly worth of being considered, i.e. the 
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parton cm. angular dependences of the differential cross section in various helicity channels. 
We have chosen four cm. energy values, y/s = 300, 500, 1000, and (academically) 2000 GeV 
and retained for sake of comparison the full angular range — 1 < cos 9 < 1 (possible angular 
cuts will be considered separately). We have only retained those terms that are numerically 
meaningful, leaving aside the invisible ones. In the Figure, for simplicity, we show only the 5 
amplitudes which are leading at high energy. These are the three asymptotic ones generated 

by the helicity amplitudes previously defined and two extra ones, corresponding to F 

and The important points to be noticed are the following ones: 

1. The relative relevance of the different helicity differential cross sections changes dras- 
tically with the scattering angle for the two lower energy points. As one sees, the 
scattering in the nearly backward region is totally dominated for y/s = 300 — 500 GeV 
by the two non asymptotic quantities; the weight of the asymptotic differential cross 
sections becomes dominant when moves to the forward direction, where the overall 
numerical size is, though, smaller than that of the backward region. 

2. Although less evidently, these features survive also at the next energy point \fs = 

1 TeV. More precisely, the size of the F distribution remains essential in the 

backward region. 

3. One might start doubting about the validity of our asymptotic assumptions. To show 
that this is not the case, we have plotted the distributions in the last sub-figure, for 
the (academic) point = 2 TeV. As one sees, the features at this energy are those 
that would expect at (sufficiently!!) high energies:the largely dominant contribution 
is that of two of the asymptotic quantities, more precisely F_ ++0 and F. 

In conclusion, we see that the contribution of the non asymptotic helicity amplitudes, for 
which no Sudakov expansion has to be expected, is essential for realistic (i.e., qualitatively 
< 1 TeV) energies. A proper asymptotic behavior seems to eventually set in, but only at 
higher energies (say, ~ 2 TeV), where the possibility of detecting the signal appears, least 
to say, debatable. Although these features were derived by an analysis performed in Born 
approximation, we expect that the complete results that will follow will be consistent with 
these preliminary impressions. 
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V. RESULTS 



The successful results of our previous tests have encouraged us to prepare with a reason- 
able amount of confidence a numerical C++ code that contains the complete tested one-loop 
expression of all the components of the considered process. This program has been called 
MINSTREL and is nowadays working and available. Thanks to this code we are now able 
to provide numerical predictions for the complete electroweak effect of the MSSM at one 
loop on the realistic observables of the associated tW~ production. 

With this aim, we have returned to Eq. ()47|) and have considered a set of SUSY bench- 
mark points that appeared to us suited for our analysis. More precisely, we have retained 
representative points whose SUSY masses values are not light ( but not even dramatically 
large) and also points whose masses are, conversely, light (in our language, lighter than, say, 
400-500 GeV). Also, we have used points whose only essential difference is the value of tan (3, 
that is allowed to become definitely large (50) in one of the two cases and still appreciable 
(10) in the second one. In this way, we should be able to compare the complete results with 
those that we found in the light SUSY Sudakov approximation. For practical reasons, we 
will only show the results of our analysis for a choice of four representative points. Two 
of them are the ATLAS DC2 SU1 and SU6 points jig ]: the remaining two are two points 



whose spectrum has been evaluated by the code SUSPECT |17[ and that we have called 
LSI, LS2 where LS stands for Light SUSY. To make the reasons of our choice evident, we 
have given in Tab. (JH) the values of the various SUSY masses, and of tan/?, that correspond 
to the four choices. One sees that the first two points correspond to a not light choice, with 
two different values of tan/?; for the last two points, a light SUSY scenario is assumed, with, 
again, two different tan (3 values. Two final technical points have to be now added: 

a) Our calculations have been performed with a value of PT,mm = 15 GeV. 

b) In the calculations, we have included a QED soft photon contribution, computed 
assuming an upper value of the soft photon energy AE = 0.1 GeV. As we anticipated, 
the full treatment of the essentially Standard Model hard photon emission will be 
contained in a dedicated paper 
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A. Unpolarized observables 



1. Effects in the distribution da/ds 

We can now show the first results of our calculations. In Fig. 0] we have drawn the 
relative effect at one loop of the MSSM, and also of the SM alone, for the four choices of 
benchmark points. The calculation stops at y/s = 1 TeV as we announced. From a glance 
at the different Figures, a number of (preliminary) conclusions can already be drawn. In 
particular: 

a) The genuine SUSY effect, i.e. the difference between the MSSM and the SM, remains 
systematically small (a relative few percent) for all choices of the benchmark points in 
the considered (realistic) energy region. In this sense, a measurement of the invariant 
mass distribution of the process does not appear to be a promising way of detecting 
genuine SUSY effects in the MSSM with mSUGRA symmetry breaking (this conclu- 
sion could be not valid for different supersymmetric models or symmetry breaking 
scenarios) . 

b) The relative effect of the considered MSSM is not, though, negligible. As one sees, it 
varies from positive to negative values in the lowest part of the region, remaining sys- 
tematically negative for larger energies and reaching a common value of approximately 
ten percent around 1 TeV. This large energy negative shift from the Born level calcu- 
lation appears a characteristic property of the considered MSSM model, independent 
of the values of the parameters that were assumed in our analysis. 

2. Ratios of partially integrated cross sections: a proposal 

Since there is a wide energy region where the one loop effects are appreciable, i.e. \/s > 
500 GeV, we can split it in two parts, compute the associated integrated cross section, and 
evaluate the ratio R of the two partial cross sections. This investigation is motivated by the 
following remarks concerning general properties of R: 

a) It should be free of several systematic experimental errors; 

b) It should be free of several QCD effects (same pdfs, same virtual corrections); 
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c) It should be essentially unaffected by photon radiation effects. 

To give an explicit numerical example, we have considered the scenario SU6, and have 
split the high energy region in two parts: 

^threshold = m t + m w < y/s < 400 GeV, > 400 GeV. (49) 

We call <t_ and a + the integrated cross section J(da/ds) ds in the two regions, and define 
R = <j + /a_. 

We denote by e± the relative MSSM effect on the two cross sections. We also denote by 
N± the expected number of events associated to the two regions. Of course, N± = £ <r±, 
where C is the luminosity. If we call Amssm-R and A stat -R the MSSM and statistical shifts 
on R we have 

A M ssm# = #(£+-£-), A statJ R = R ( + ) . (50) 



'+ 

In our test case, the Born value is R ~ 0.58 and the difference e + — e_ gives a shift 0.57 — > 0.60 
of about 3.5 %. the purely statistical error computed with a luminosity C = 10 ftT 1 gives 
A s tat-R — 0.002, i.e. a shift about 10 times smaller than the MSSM effect. 

We conclude that radiative effects in ratios like R are beyond the statistical noise. Of 
course, systematic errors are expected to dominate over statistical ones. Thus, a detailed 
dedicated experimental study of the process reconstruction will be crucial to assess R as a 
realistic observable and a potential precision test of the electroweak sector of the considered 
MSSM. 



3. Sudakov-like parameterizations 

To conclude the unpolarized session, we have tried to give an effective parametrization of 
the full one loop effect in the spirit of the logarithmic Sudakov expansion. As we remarked, 
a straightforward comparison with the results described in Sec. Illll is hampered by a variety 
of problems, that we now emphasize: 

a) Box diagrams are functions of the Mandelstam invariants t, u beside s. At small or 
large angles, these can be small (compared to the internal squared masses and s) and 
spoil the validity of the Sudakov approximation. 
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b) At high but moderate energies (below 1 TeV) there are several subleading helicity 
channels which are relevant and non negligible. These channels certainly admit a 
Sudakov expansion. However this is not as simple as that of the leading channels. The 
coefficients of the expansion for these amplitudes have not been investigated before 
and could or could not turn out to be simple combination of quantum numbers and 
couplings as in the leading case. 

c) In the MSSM, we have sparticles with masses around 300-400 GeV, even in the lightest 
considered scenario LS2. The extent to which they can be regarded as small can only 
be determined by an explicit numerical comparison of the two calculations. Indeed, by 
a careful inspection of the various involved diagrams one sees that box diagrams can 
display a rather delayed asymptotic behavior. In practice, if the typical virtual masses 
are of order m, there are box diagrams with asymptotic behavior ~ log(y/s/m') where 
m! can be 4-5 times larger than m, depending in particular on the scattering angle. 
This large effective scale contributes a large energy independent constant shift in the 
difference between the Sudakov and the one loop calculations. Also, since we always 
require \fs 3> m', it pushes forward the energy range where the expansion is accurate. 

As a consequence of these remarks, the difference between the full one-loop MSSM effect 
and the Sudakov approximation is expected to be a small, slowly varying function of the 
energy, at least in the considered energy range. On the contrary, in the Standard Model, 
all masses are quite light compared to the typical 500-1000 GeV parton energy and we can 
hope to observe a better accuracy of the Sudakov expansion. 

All these expectations are confirmed by actual calculations. As an illustration, we show in 
Fig. (0) the comparison between the full one loop and the Sudakov calculations of the effect 
in the distribution da/ds. The left panel shows the Standard Model case. The right panel 
shows the LS2 MSSM scenario, which is the lightest considered. For purpose of comparison, 
we have switched off QED radiation and set M 7 = Mz- We have computed the effects up to 
unrealistic values (2 TeV) of the energy, just to emphasize the convergence at high energy. 
The Sudakov approximation is evaluated with a common scale M in the double and single 
logarithms. The best value of M is an important issue and will be discussed below. The 
main features of the figure are the following. 

a) In the Standard Model, we choose M = Myy. We observe a remarkable agreement. 
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The expansion is rather accurate down to energies y 7 ? ~ SOOGeV^. The relevant 
scale is the electroweak breaking one ~ M\y and there are no large constant (energy 
independent) contributions. 

b) In the LS2-MSSM case, we adjust M in order to have the same slope in the two curves. 
We have found the optimal value M = 120 GeV. With this choice, there is a large 
but constant shift of about +6% with respect to the Sudakov calculation as shown 
in the upper right panel where we show the two curves. To emphasize the energy 
independence of the shift, we show in the lower right panel the same one loop curve 
together with a shifted version of the Sudakov one which has been moved upward by 
a constant +6%. The agreement is again remarkable, exactly like in the Standard 
Model case. 

In principle, these features could be useful if one were interested in preparing a complete 
NLO parametrization of the process, that includes QCD effects and decay simulation by 
Monte Carlo. The expected smoothness of the radiative effects beyond thresholds can be 
exploited to replace the full calculation by simple (model dependent) interpolating expres- 
sions. This is particularly relevant in the SM case where the Sudakov-like parametrization 
of the process is fixed and does not depend on any model parameter, but only on the kine- 
matical cuts. 

This remark concludes our presentation of the unpolarized effects. We move now to the 
discussion of the possible polarized observables of the process. 

B. Polarized observables 

1. Final top asymmetry: one loop effects in A^r 

A special property of the tW~ production process is the fact that, in principle, the 
polarization of the final top quark and/or W boson can be measured. This fact, that was 
first considered in a previous reference [la], leads to the introduction of new observables, that 
we shall try to list and to discuss in what follows. The first possibility is that of measuring 
the final top polarization. In the process that we are considering, the final top can have in 
principle both helicities, as one can see from the expressions of the helicity amplitudes given 
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in Section 3. In correspondence to the two possible choices, we shall define two different 
differential cross sections, that we shall define as daL,R/ds, that are the analogues of Eq. (|37jl 
where only the contributions from the two types of final top have been retained. Plotting 
these quantities at variable y/s, as we did for the total unpolarized cross section, would 
lead to conclusions that do not much differ from those already given in the previous part 
of this Section: the genuine SUSY effect is still rather modest. Again, the overall MSSM 
effect is, though, not small. This could be seen in the plots of the two distributions, but 
from our previous discussion we believe that it might be preferable to consider, again, ratios 
of cross sections. With this aim, we have defined the ratio of the integrated cross sections 
asymmetries, i.e. the quantity 



A L r{s) = — r -— i rT , with a L>R {s) = / — — '— ds'. (51) 

<W) + oji(s) J E 2 ds' 

threshold 



Figures (|6I7|) shows the values of Alr at variable \/s. One sees that, considering a realistic 
value e.g. yfs = 500 GeV, the one-loop effect on the asymmetry reaches in all considered 
SUSY scenarios an absolute value of slightly less than 1 %. This number should be compared 
to the realistic overall uncertainty. For the reasons that we have discussed previously, we 
expect essentially a dominance of the purely statistical experimental error, whose size will 
depend on the available integrated luminosity. Lacking a dedicated experimental analysis 



(in preparation), we can use as a guidance the pre 
single top production process, the t-channel one) of 
percent. 



iminary quoted value (for a different 
that is a (mainly statistical) four 



2. Final W asymmetry: one loop effects in Atl 

In the tW~ production process, the final W~ is real. Therefore one can, in principle, 

measure the primary W polarization. Assuming that this is the case, we have defined two 

quantities that are the analogues of Eq. (J5Tj) and, starting from them we have introduced 

the transverse-longitudinal asymmetry, defined as 

a Wr (s) - o- Wl (s) f s da WT 

A TL (s) = -— — , with a w (s) = / - - ds . 52 

o- Wt (s) + o- Wl {s) J E 2 ds' 



threshold 



The numerical values of Atl are shown in Figs. (|8|9|) . In all cases the one-loop effect at the 
point yfs = 500 GeV has an absolute value of about 0.5 %. We do not have yet at disposal 
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a suitable experimental analysis for this asymmetry, that is in fact being carried on 1 



VI. CONCLUSIONS 

In this paper, we have performed the first complete electroweak one-loop analysis of 
the associated tW production process in the MSSM with mSUGRA mechanism of SUSY 
symmetry breaking. This has been done using a numerical program, MINSTREL, that 
satisfies the three constraints of cancellation of ultraviolet and infrared divergences and of 
reproduction of asymptotic Sudakov expansions. We have considered various experimental 
potential observables, both for unpolarized and for polarized production. We have found a 
relatively small genuine SUSY effect for the representative SUSY benchmark points that we 
have selected, and a possibly appreciable, mostly of SM origin, overall one- loop effect. We 
have proposed a number of new observables, in general ratios of experimentally measurable 
quantities, that would be essentially free of disturbing theoretical QCD and experimental 
systematic uncertainties. For these quantities, the predictions of the MSSM would be rather 
precise, making them appear as possible precision tests of the involved genuine electroweak 
content of the model. The extension of our results to a different MSSM scenario or to differ- 
ent SUSY models would be straightforward. The still missing corresponding experimental 
analysis of the various proposed observables of the process is being carried on, and will 
appear in a more complete forthcoming paper 
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TABLE I: Table of spectra for the various benchmark points. All entries with the dimension of a 
mass are expressed in GeV. The spectra have been computed with the code SUSPECT Qj. 
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FIG. 1: Born diagrams for the process bg — > iVF . 
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FIG. 2: Generic diagrams for the process bg — > tW . We list only the vertex corrections and the 
box diagrams. The labels S, F, and V denote generic particles with spin 0, 1/2, and 1. 
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FIG. 3: We show the energy and angular dependence of the 5 helicity amplitudes which are leading 
at high energy. Of course, these include the three amplitudes which are not mass suppressed. In 

addition, we show the next relevant amplitudes which are the mass suppressed ones F o an d 

An inspection of the Figure shows that below 1 TeV the mass suppression is not effective, 

especially in the backward region. 
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FIG. 5: Comparison between the full one loop and Sudakov calculations of the effect in the 
distribution da/ds. The left panel shows the Standard Model case. The right panel shows the 
LS2 MSSM scenario, which is the lightest considered. In the upper figure we show the two curves, 
whereas in the lower panel, we have shifted the Sudakov effect by a constant +6%. For purposes 
of comparison, we have switched off QED radiation and set M 7 = Mz- 
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